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ABSTRACT 


Based on gas-phase laboratory spectra at 6K, Campbell et al. (2015) con 


firmed that the diffuse interstellar bands (DIBs) at 9632.7 and 9577.5 A are due 
to absorption by the fullerene ion CgQ- They also reported the detection of two 
other, weaker bands at 9428.5 and 9365.9 A. These he in spectral regions heavily 
contaminated by telluric water vapour lines. We acquired CFHT ESPaDOnS 
spectra of HD 183143 close to the zenith and chopped with a nearby standard 
to correct for the telluric line absorption which enabled us to detect a DIB at 
9365.9 A of relative width and strength comparable to the laboratory absorption. 
There is a DIB of similar strength and FWHM at 9362.5 A. A stellar emission 
feature at 9429 A prevented detection of the 9428.5 A band. However, a CFHT 
archival spectrum of HD 169454, where emission is absent at 9429 A, clearly 
shows the 9428.5 A DIB with the expected strength and width. These results 


further conhrm C^q as a DIB carrier. 


Subject headings: ISM: general — ISM: lines and bands — ISM: molecules 
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Introduction 


Hundreds of diffuse interstellar bands (DIBs) are seen in the spectra of stars dimmed 


and reddened by interstellar dust (Herbig 1995). There have been many unsuccessful 


attempts to assign neutral or charged molecular carriers to either individual or families of 


these bands. Foing and Ehrenfreund (1994) proposed that a pair of DIBs which they had 


detected at 9632 and 9577 A were due to the fullerene ion CgQ because of their proximity to 
laboratory absorption wavelengths observed in a neon matrix (Fulara et al.||1993). Campbell 


et ah (2015) have now proved that these DIBs are indeed due to CgQ. This is based on 


measurement of not only the wavelengths, but also the FWHM and relative intensities of 
these bands in the gas phase at 6 K. This successful identihcation, which came long after 
Kroto first predicted the presence of C^q in the interstellar medium ( |Kroto 1987), may 
be regarded as a first step towards unlocking the 90 year old mystery of the DIB carriers 


(Herbig (1995), Snow and McCall|| (2006)). 


Campbell et ah (2015) also reported two weaker bands of roughly equal strength (3:2) 


at 9428.5 and 9365.9 A in their laboratory spectra (their Table 1). Unfortunately, for ground 
based observations, both of these wavelengths lie in spectral regions heavily contaminated 
by telluric water vapour (WV) lines as can be seen in Figure The standard technique for 
their removal is division of the reddened star spectrum by one from an unreddened star of 
similar spectral type observed through the same air mass. 

In this paper, we report observations of HD 183143 chopped with those of a similar 
standard lying within 2° on the sky. Atmospheric WV absorption is likely to change 
with both time and direction as well as with air mass which is the reason for frequent 
chopping. HD 183143 is a bright, reddened, B7 supergiant which exhibits all of the known 


DIBs with considerable strength (Herbig 1995) including those at 9632 and 9577A. It also 


passes through the zenith at CFHT where the atmospheric precipitable water level can be 
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significantly lower than for sites at lower altitnde. 


2. Observations of HD 183143 


We were granted two hours of CFHT Director’s Discretionary Time with ESPaDonS 


(Donati 2003) on the night of 2015-07-28 UT. HD 183143 and a nearby standard of similar 
spectral type, HR 7437, were observed within one hour of meridian crossing with a sky 
temperature of ~ —38 °C indicating a low water vapour level. ESPaDonS was configured 
in spectropolarimetric mode consequently, a single spectrum consisted of four individual 
sub-exposures taken at different polarimeter retarder configurations. All of the spectra were 


reduced at CFHT with the Upena pipeline using the Libre-ESpRIT code (Donati et al. 


1997). The spectral resolution is approximately 65,000 or 0.05 A per pixel. Details of the 


stars are given in Table 1 which includes V and I magnitudes, and colour excess, E(^b-v)- 


Observations alternated between the two targets. In this way, four sets of spectra were 
acquired with air masses differing by < 0.01 within each pair. The total exposure times of 
600 and 720 s within each set for HD 183143 and HR 7437, respectively, were chosen to give 
similar signal levels for both stars in the 9400 A region. 

Telluric lines in typically 20 A regions of interest in the spectrum of HD 183143 
were removed using the subsequent HR 7437 spectrum as standard in the IRAhj^ telluric 
program. The four corrected spectra of HD 183143 from each of these spectral regions were 
then combined to produce the single high S/N spectra in Figure]^ Estimation of S/N is 
complicated by the presence of the telluric lines within which S/N is low. Based on our 


^IRAF is distributed by the National Optical Astronomy Observatory, which is oper¬ 
ated by the Association of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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calculations, S/N ~ 1500 A ^ (~ 330 per pixel) for HD 183143 in the region of the 9577.5 A 
line after WV removal. 

There is a short gap in the ESPaDOnS spectral coverage between 9608 and 9636 A 
where the order lies off the detector. This gap prevented our looking for the 9632.7 A band. 


3. Results 

The wavelength region of interest from 9300 to 9600 A is covered in a single echelle 
spectral order (24) and shown in the lower plot of Figure for HD 183143. All of the 
strong lines, some near saturation, are due to telluric WV. The upper plot is the result 
of division by the spectrum of HR 7437. Note that the telluric line cancellation has not 
been optimized in the hgure since in this case a single correction was applied to the full 
300 A region displayed. The narrow stellar hydrogen P8 line at 9546 A in HD 183143 is 
superimposed on the reflex of the much broader P8 line in HR 7437 illustrating the impact 
of absorption lines in the standard star. The strong 9577.5 A band of CgQ is easily seen. 
Several other DIBs are also visible. 

To optimize telluric line removal we performed these corrections over more limited 
spectral ranges containing the wavelengths of the C^q laboratory absorptions. This produced 
higher S/N spectra for the 9577.5, 9428.5 and 9365.8 A regions as well as a region centered 
on the interstellar potassium line at 7699.0 A in order 29. 

Small zero point adjustments were made to the levels of the HD 183143 and HR 7437 
spectra to minimize residuals in the ratios at the centres of the strongest WV lines. The 
ratios are shown for the four regions in Figure The persistence of some residuals at the 
centres of the strongest telluric lines is likely associated with non-linearity of the analog to 
digital converter over the extreme range of intensity of the lines rather than inequality of 
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water vapour optical depth. 

Figure 1^ also shows the two components of the HD 183143 interstellar 7699.0 A 
potassium absorption line after removal of nearby telluric lines. The stronger component 
has a radial velocity of —12 km s“^, the other component +4 km s“^. At 9577.5, 9428.5 
and 9365.8 A, —12kms“^ corresponds to a displacement of —0.4 A. We use this value in 
what follows assuming that the corresponding cloud is the main contributor to the detected 
DIBs. 


Campbell et ah (2015) fitted Gaussians to the laboratory absorptions at 6K and 


listed wavelengths and FWHMs (their Table 1). The DIBs in Figure]^ were also fitted 
with unconstrained Gaussians as shown by the dashed lines including the obvious DIB at 
9362.5 A. No DIBs could be observed in the region of the 9428.5 A laboratory absorption 
because of the emission feature at 9429 A. The emission must be associated with HD 183143 
because the comparison star, HR 7437, is a rapid rotator (usini = 265kms“^) such that 
absorption lines in HR 7437 would produce much broader (> 8A) apparent emission features 
in the ratio. There are many such emission lines between 9300 and 9600 A for HD 183143. 


Jenniskens et ah (1997) published spectra in this region, and encountered emission lines in 


either their reddened or standard stars (their Figure 4). It is interesting to note that in the 
latter study, which was published long before the gas phase spectrum of G^g recorded, 
both a ‘depression’ near 9428 A and the suggestion of a weak feature around 9366 A were 
reported. 


Figure shows the regions of the 9577.5, 9428.5 and 9365.9 A laboratory bands for the 
B1 la star HD 169454 divided by the spectrum of the rapidly rotating (usini = 317kms“^) 
AO V star HD 177724, again using the IRAF telluric program. Details of the stars are 
given in Table 1. These processed spectra were downloaded from the Ganadian Astronomy 
Data Gentre’s GFHT archive. Both stars were observed in the “object only” spectroscopic 
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configuration of ESPaDonS on 2005-05-21 UT. We selected HD 169454 because it has 
none of the emission features found in HD 183143. Unfortunately, the standard was not 
observed at an identical air mass and so there are large WV line residuals at 9427.5 and 
9428.2 A which have been omitted and replaced by averages in the plots and when htting a 
Gaussian. The principal interstellar K line (7699.0 A) is unresolved with a radial velocity 
of —11 km s“^ which translates to —0.4 A at 9400 A. Table 2 summarises the results of 
unconstrained Gaussian hts to the features in all of the ratioed spectra. 

Given the weakness (~ 1%) of the detected bands and the presence of apparent 
emission features, setting of continua was subjective, particularly in the case of HD 169454. 
Further, the GGD detector is partially transparent at such long wavelengths which generates 
interference fringes. These are not well cancelled in the flat helding process but the 
additional normalization when dividing by the telluric standard considerably reduces their 
amplitude. Gonsequently, it is hard to establish values for systematic errors with conhdence 
particularly for the FWHM in the case of HD 169454. In the table, only formal htting 
errors are given. 


4. Comparison with CgQ laboratory bands 

The agreement of the astronomical and laboratory wavelengths is remarkably good 
considering the somewhat ad hoc nature of assigning an interstellar velocity. In addition, 
support for the assignment to G^q comes from comparison with the relative intensities 
reported in the laboratory measurements. The depth of the astronomical bands (Table 
2) can be used to estimate the relative intensity of the absorption bands for symmetric 
prohles. The observations towards HD 183143 give a ratio of 1:0.3 for the 9577.5 and 
9365.9 A bands. This matches with the laboratory ratio of 1:0.2 which was reported with 


an estimated uncertainty of around 20 % (Gampbell et al. 


2015). Furthermore, the ratio 





of the depth values derived from the archive HD 169454 spectra, 1:0.3:0.2 for the 9577.5, 
9428.5 and 9365.9 A DIBs, are also in agreement with 1:0.3:0.2 measured in the laboratory 


(Campbell et al. 2015). 


The FWHM of 3.3 A for the interstellar 9577.5 A band in HD 183143 is comparable to 


the value of 2.85 A reported by Foing and Ehrenfreund (1997) and 3.0 A by Jenniskens et 


al. (1997) for other reddened stars. In the laboratory study at 6K the measured FWHM 


are 2.5 A. At this low temperature, the FWHM of the C^q rotational prohle is around lA 


(Edwards and Leach 1993). This indicates that the laboratory bands are broadened by the 
2ps lifetime of the excited electronic state. Only at temperatures above 30 K would one 
see a further broadening of the rotational profile. The FWHM values for HD 169454 are 
not robust because of incomplete correction for telluric contamination and arbitrariness in 
setting the continuum. A more likely source for the extra broadening of the 9577.5 A band 
in HD 183143 would be near coincidence with another DIB. For example, there is a close 
DIB to the 9365.9 A band at 9362.5 A (Figure [^. 


5. Astromomical implications 

The conhrmation of the presence of C(|q in diffuse clouds raises the question of its 
interstellar abundance. This has been considered previously, see, for example ( Herbig||2000 ), 
and references therein. The calculation of column density relies on the knowledge of the 
individual band oscillator strength, /(A), the wavelength of the absorption band maximum, 
Ac, and the equivalent width, W. For this purpose the value for the electronic transition 


which was estimated in the matrix isolation spectroscopy study by Fulara et al. (1993), 


fe = 0.003 — 0.006, has invariably been used. However, this value is an order of magnitude 


smaller than indicated by theory (Bendale et al. 1992). We believe that the value from the 


matrix study is underestimated and in the following /g = 0.05 is used. 
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In order to evaluate the /(A) value for the 9577.5 A band the electronic oscillator 
strength must be weighted by the Franck-Condon factor. In the case of C^q, most of 
the oscillator strength is localized in the two strong features at 9632.7 A and 9577.5 A 


(Campbell et ah 2015). Based on observation of only weak CgQ laboratory absorptions to 
shorter wavelength than 9428.5 A, we estimate the 9577.5 A band to account for ~ 30 % 
of the intensity and therefore obtain /(9577.5A) ~ 0 . 02 . For the 9577.5 A DIB observed 
towards HD 183143 the equivalent width is estimated to be W 0.3 A (Table 2 ). Using 
these values gives a column density, A^(C^q), of around 2 x lO^^cm”^. For comparison, this 
is an order of magnitude smaller than that of H;^ but similar to CH”^, both of which were 


observed towards HD 183143 (McCall et ah 2002). 


6. Conclusion and outlook 


Two weak diffuse interstellar bands have been detected with wavelengths, FWHM 


and relative intensities in agreement with the laboratory spectrum of CgQ (Campbell et ah 


2015). These results provide further compelling evidence for the presence of this molecular 


ion 


in the interstellar medium. It has been suggested that the 9632.7 A and 9577.5 A 


absorptions, and hence CgQ, are also present in protoplanetary nebulae (Iglesias-Groth and 


Esposito 2013). The existence of the neutral fullerenes Ceo and C 70 in a young planetary 


nebula (Cami et ah 2010) and reflection nebulae (Sellgren et al. 2010) has been conhrmed 


through their infrared (IR) transitions. There is also a report of IR emission features from 


a reflection nebula which have been proposed to be due to CgQ (Berne et al. 2013). The 


identihcation of CgQ in diffuse clouds leads to intriguing questions regarding the role of 
fullerenes with respect to the formation of smaller carbon based molecules such as those 
identihed in dense interstellar clouds by radioastronomy. 























Table 1. CFHT ESPaDonS Observations 


star 

Sp/L 

1 

V 

E(b-v) 

HD183143 

B7 la 

4.79 

6.86 

1.28 

HR 7437 

B8 Illn 

5.10 

5.00 

0.00 

HD 169454 

HI la 

5.13 

6.71 

1.12 

HD 177724 

A4 IV-Vnn 

2.98 

2.99 

0.00 


Table 2. Laboratory and Interstellar C^q Bands^ 


laboratory 

Ac f^rel^ 

FWHM A 

Ac* A 

HD 183143 

depth % 

FWHM A 

Ac* A 

HD 169454 

depth % 

FWHM A 

9577.5 ±0.1 

1 

2.5 ±0.2 

9577.4 ±0.02 

9.1 ±0.05 

3.3 ±0.04 

9577.2 ±0.03 

4.0 ±0.06 

3.5 ±0.06 

9428.5 ±0.1 

0.3 

2.4 ±0.2 

- 

- 

- 

9428.4 ±0.1 

1.2 ±0.05 

3.2 ±0.1 

9365.9 ±0.1 

0.2 

2.4 ±0.2 

9365.7 ±0.02 

2.4 ±0.03 

2.5 ±0.04 

9365.6 ±0.1 

0.7 ±0.1 

2.1 ±0.2 


^central wavelengths (Ac), depths, FWHM and errors from Gaussian fits to the lab and astronomical bands. No 
systematic errors included. 

1 relative cross section 

^corrected by +0.4 A for the interstellar K line off-set 
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I O' _!_^_ ni I n _ I ' I ' _^^_n_j^_M__!_J__| 

9300 9350 9400 9450 9500 9550 9600 


Wavelength (A) 

Fig. 1.— The lower spectrum of HD 183143 taken at CFHT contains the region of the C^q 
absorption bands observed in the laboratory. The spectrum is dominated by strong telluric 
water vapour absorption lines. The upper plot follows division of the lower HD 183143 
spectrum by that of HR 7437 which eliminates most of the water vapour lines. The feature 
at 9546 A is caused by hydrogen (P8), the narrow absorption component is from HD 183143 
while the broad reflex (apparent emission) is from the broader P8 absorption line in HR 
7437. Weak emission features from HD 183143 can be seen throughout the upper plot as 
well as a number of DIBs. The C)|q 9577.5 A band is clearly visible. 
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Wavelength (A) Wavelength (A) 


Fig. 2.— The spectral regions of interest for HD 183143 after cancelation of tellnric water 
vapour lines by division with HR 7437 and adjustment of zero levels. The vertical dotted 
lines indicate laboratory rest wavelengths for the 7699.0 A potassium line, and the absorption 
bands of C^q. The dashed lines are Gaussian hts for which parameters are listed in Table 2. 
The emission feature at 9429 A prevented any interstellar band detection for HD 183143. A 
DIB is also seen at 9362.5 A and fitted separately with a Gaussian. 
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Wavelength (A) Wavelength (A) 

Fig. 3.— The same spectral regions for HD 169454 after cancelation of telluric water vapour 
lines by division with HD 177724 and adjustment of zero levels. The vertical dotted lines 
indicate laboratory rest wavelengths for the 7699.0 A potassium line, and the absorption 
bands of Cgg- The dashed lines are Gaussian hts for which parameters are listed in Table 2. 
In the panel for 9428.5 A, large water vapour line residuals at 9427.5 and 9428.2 A have 
been omitted and replaced by averages. The DIB seen at 9362.5 A is htted with a separate 


Gaussian. 

















